Trace formula for chaotic dielectric resonators tested with microwave experiments 
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We measured the resonance spectra of two stadium-shaped dielectric microwave resonators and 
tested a semiclassical trace formula for chaotic dielectric resonators proposed by Bogomolny et al. 
[Phys. Rev. E 78, 056202 (2008)]. We found good qualitative agreement between the experimental 
data and the predictions of the trace formula. Deviations could be attributed to missing resonan- 
ces in the measured spectra in accordance with previous experiments [Phys. Rev. E 81, 066215 
(2010)]. The investigation of the numerical length spectrum showed good qualitative and reasonable 
quantitative agreement with the trace formula. It demonstrated, however, the need for higher-order 
corrections of the trace formula. The application of a curvature correction to the Fresnel reflection 
coefficients entering the trace formula yielded better agreement, but deviations remained, indicating 
the necessity of further investigations. 

PACS numbers: 05.45.Mt, 42.55.Sa, 03.65.Sq 



I. INTRODUCTION 

The development of dielectric microresonators with 
their manifold applications [IH2] like microlasers [4] 
has triggered much interest in their theoretical descrip- 
tion. Dielectric resonators are governed by the vecto- 
rial Helmholtz equation, whose exact treatment is gener- 
ally not achievable. However, semiclassical methods have 
been applied successfully to the modeling of, e.g., the 
far-field patterns of microlasers [SHI] or the localization 
of their resonant modes, so-called scarring [8HT5]. These 
involve the periodic orbits (POs) of the corresponding 
classical system, a two-dimensional (2D) dielectric bil- 
liard in the case considered here. The POs are, for exam- 
ple, connected to the density of states via trace formulas 
[l4Tfl"6] . and thus, also to the spectral properties of the 
corresponding resonator. Trace formulas for 2D dielec- 
tric resonators have been proposed [TTJ [T8] and tested 
experimentally with microlasers |19L I20j and microwave 
resonators |21H23j for various geometries, which predomi- 
nantly correspond to dielectric billiards with regular clas- 
sical dynamics. 

The aim of the work presented here has been a thor- 
ough experimental test of the trace formula for chaotic 
2D dielectric resonators. We used two microwave res- 
onators in the shape of stadia with different aspect ra- 
tio. The stadium billiard is a fully chaotic system [24] 
and has been investigated thoroughly theoretically j25j- 
134] and experimentally 35 42] in the context of quantum 
chaos. Stadium-shaped microlasers have been studied ex- 
perimentally in Refs. p3Hl3l USS] , and in Ref. [20], the 
length spectrum of a stadium-shaped microlaser was de- 
termined from measurements. However, there still was 
lack of a detailed comparison of the experimental length 



spectrum of a chaotic dielectric cavity with the trace for- 
mula prediction of Ref. [17 . This has been the motiva- 
tion of the present work. 

The article is organized as follows. Section [H] describes 
the experimental setup and the measured frequency spec- 
trum, Sec. |TlT] summarizes the salient features of the trace 
formula for chaotic dielectric resonators, and Sec. HVI con- 
tains the comparison of the experimental length spectra 
with the trace formula predictions. In Sec. [V] the mea- 
sured data and the trace formula predictions are com- 
pared to numerically calculated data, and a correction to 
the trace formula at curved boundaries is investigated. 
Section [Vl] closes with some concluding remarks. 



EXPERIMENTAL SETUP AND 
FREQUENCY SPECTRUM 



II. 



Two flat Teflon plates in the shape of a Bunimovich 
stadium were used as microwave resonators. The ge- 
ometry of the plates is shown in Fig. |l(a)| The radius 
of the semicircle of the first stadium, designated SI, is 
Ri = 150.1mm and the length of its straight edge was 
L\ = 200.2 mm. Its aspect ratio is thus ei = Li/(2i?i) = 
0.67. It had an index of refraction of ri\ = 1.425. The 
second stadium, S2, had the parameters R 2 = 149.9 mm, 
L 2 = 99.8 mm, e = 0.33 and n 2 = 1.404. The index of 
refraction was determined in both cases with the same 
technique as in Ref. [5T]. The small difference between 
Ri and R 2 is due to manufacturing uncertainties. The 
frequency range of interest, 5-20 GHz, corresponds to 
kR = 15.7-62.9, where k is the wave number. Both 
Teflon plates had a thickness of d = 5 mm. The dielectric 
plates were put between two copper plates as shown in 



Fig. 1 1(b) Below the frequency [4"T] 
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FIG. 1. [(a)] Geometry of the two Teflon plates SI (left) and S2 
(right). The solid lines indicate the boundaries of the plates 
and the dashed lines the common border of the semicircular 
and the rectangular parts of the stadia, where R is the radius 



of the semicircles and L the length of the rectangular part, (b) 
Schematic side view (section) of the experimental setup (not 
to scale) . The Teflon plate with thickness d is put between two 
copper plates. Two antennas are lead through small holes in 
the top copper plate next to the sidewalls of the Teflon plate. 
Reprinted from Ref. [21] . 



only two-dimensional transverse magnetic modes exist in 
the resonator, where c is the speed of light in vacuum 
and n is the index of refraction. In the cases considered 
here, feu ~ 20 GHz. The setup is similar to that used in 
Ref. [21 j . Two vertical wire antennas were used to cou- 
ple microwave power into and out of the resonator. They 
were placed next to the sidewalls of the Teflon plates as 
shown in Fig. |l(b) and protruded 1 mm into the space 
between the copper plates. A vectorial network analyzer 
(PNA 5230A by Agilent Technologies) was connected to 
the antennas via coaxial rf cables. It measured the com- 
plex transmission amplitude Sb a {f) from antenna a to 
antenna b for a given frequency /. Its modulus squared 
is given by 



\Sba\ 



P, 



b. out 



Pa- 



(2) 



where P a ^ n is the power coupled into the resonator via 
antenna a and Pb,ou% that coupled out via antenna b. The 
plot of \Sb a {f)\ 2 versus the frequency yields the transmis- 
sion spectrum of the resonator. 

Figure [2] shows the measured frequency spectrum of sta- 
dium S2. It comprises a multitude of resonances of di- 
verse widths on top of a slowly oscillating background 
that is attributed to direct transmission processes be- 
tween the antennas. The quality factors Qj = fj/Tj of 
the resonances, where fj is the resonance frequency and 
Tj the full width at half maximum of resonance j, are in 
the range of 100-1200, with a mean value of (Q) ps 300. 
The relatively low quality factors are attributed to the 
large radiation losses of the stadia. The resonances are 



additionally broadened due to absorption in the Teflon, 
Ohmic losses in the copper plates and the power cou- 
pled out by the antennas. The resonance frequencies 
and widths were determined by fitting Lorcntzians to the 
measured spectrum. The frequency spectrum includes a 
series of apparently equidistant resonances that are indi- 
cated by the arrows in the bottom panel of Fig. [2] Similar 
structures were observed in the lasing spectra of stadium- 
shaped polymer (n = 1.5) microlasers [201 EH E3 ■ The 
free spectral range (FSR) of this series, i.e., the frequency 
spacing between adjacent modes, is plotted in Fig. [3j The 
FSR is indeed constant over the whole frequency range 
considered here with a mean value of (FSR) = 194.6 MHz 
within a standard deviation of AFSR = 3.7 MHz. The 
fluctuations of the measured FSR can at least partly be 
attributed to the error in the determination of the re- 
sonance frequencies when modes strongly overlap, as is 
generally the case for the spectra presented in this work 
(cf. Fig. [2]). Other, systematic reasons for the fluctua- 
tions of the FSR are not known, but cannot be excluded. 
Also for stadium SI such a series of equidistant reso- 
nances was found with (FSR) = (165.4 ±3.8) MHz. This 
leads to the presumption that the corresponding resonant 
modes are localized around one or several POs (so-called 
scarring [48]), as has been observed for various dielectric 
cavities [8HT3"]. The length of the underlying PO(s) can 
be determined from the FSR via [T9"] 



/ 
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This approach is commonly used to identify the POs that 
play the dominant role for the modes of a microlaser, 
though a clear identification is not always possible [U |4"9T - 
From Eq. ^ we obtain l x = (1.272 ± 0.029) m for 
stadium SI and l 2 = (1.097±0.021) m for stadium S2. We 
will call I the path length in the following and compare 
it to the lengths of the various POs in Sec. IIVI The 
path length li corresponds to 94.8% of the circumference 
of stadium SI and l 2 to 96.1% of that of stadium S2, 
indicating that the modes belonging to the equidistant 
series are of the whispering gallery type. This supposition 
is supported by the fact that the associated resonances 
are the ones with the highest quality factors. 



III. HELMHOLTZ EQUATION AND TRACE 
FORMULA 

Flat microlasers or dielectric microwave resonators can 
only be treated approximately as 2D systems by intro- 
ducing an effective index of refraction. This approxi- 
mation, however, has only a limited accuracy as dis- 
cussed in Refs. [531 IM]- Therefore, the Teflon plates 
were squeezed between two metallic plates as depicted 
l(b)| Such resonators are described by the 2D 
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scalar Helmholtz equation [47] 
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FIG. 2. The measured frequency spectrum of the stadium S2. The inset in the top part shows the shape of the stadium with 
the positions of the antennas indicated by crosses. The arrows in the bottom panel indicate a subset of resonances that are 
approximately equidistant, i.e., they have a constant free spectral range (FSR). 
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FIG. 3. Measured free spectral range (FSR) of the subset 
of equidistant resonances (marked by arrows in the bottom 
panel of Fig. [2| of resonator S2 with respect to the frequency. 
The horizontal line is the mean value of the measured FSR, 
and the gray bar indicates the standard deviation. 



an increase of Tj . 

Trace formulas relate the density of states (DOS) of a 
wave-dynamical system to the POs of the corresponding 
classical system [TUfTS] . For an open dielectric resonator 
the DOS is 



Im (kj) 



Re (kj 



[Im(%)F 



(5) 



It can be decomposed into a smooth, average part p(k), 
and a fluctuating part pa U c(k), i-e., p — p + pa,uc- The 
smooth part is related to the area A and the circum- 
ference U of the resonator via the Weyl law given in 
Ref. [17]. In the semiclassical limit, the fluctuating part 
Pfiuc can be written as a sum over contributions from the 
POs of the corresponding dielectric billiard. For a fully 
chaotic billiard like the stadium, i.e., a billiard for which 
all POs are unstable and isolated, the trace formula pro- 
posed in Ref. |17) is based on the well-known Gutzwiller 
trace formula [III US]. It is 
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for frequencies / < f2D- Here, n(r) is the index of 
refraction at position r, k = 2nf /c is the wave num- 
ber, and E z is the z component of the electric field 
strength. Outgoing wave boundary conditions are im- 
posed on E z to account for the openness of the resonator 
PP. Therefore, the eigenvalues kj of Eq. Q are com- 
plex, where fj = cRe (kj) /(2ir) is the frequency and 
Tj = — clm(fcj) /n is the width of the resonance j (see 
Fig. [2| . The width Tj accounts for the radiation losses of 
the cavity. In a microwave experiment, the power coupled 
out of the resonator by the antennas and that absorbed 
by the dielectric material and the metal plates results in 



with the contribution of the rth repetition of the primi- 
tive periodic orbit p being 



Pp,r(k) 
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Here, £ p is the length of the primitive periodic orbit, 
M p is the monodromy matrix characterizing its stabil- 
ity, R p = flj Tj is the product of the Fresnel reflection 
coefficients Vj for each reflection of the primitive PO at 
the boundary, and p p is the Maslov index counting the 
number of focal points and caustics. The essential differ- 
ence to the Gutzwiller trace formula is an additional fac- 
tor R p accounting for the dielectric boundary conditions. 



4 



Since the polarization of the electric field, E = E z e z , is 
perpendicular to the plane of incidence of the POs, the 
Fresnel reflection coefficients are given by 55J 



n cos Oj — * / 1 — n 2 sin 2 9j 
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where 6j is the angle of incidence with respect to the 
boundary normal for the jth reflection. The POs of the 
billiards were obtained with the algorithm of Ref. [56] . 
their monodromy matrices M p were calculated following 
Ref. [57] , and the Maslov indices p p using Ref. [55] ■ The 
trace formula [Eq. |6])] is only valid for isolated, unstable 
POs. In the stadium billiard, a family of marginally sta- 
ble POs which bounce between the parallel parts of the 
billiard (bouncing ball orbits, BBOs) exists that needs 
a separate treatment [3DHS3ES]- In our case, however, 
these orbits are of little relevance since the reflection co- 
efficients for vertical incidence are quite small. 
Instead of studying the DOS itself, the Fourier transform 
(FT) of its fluctuating part was considered, 



dkp Ruc (k)e 



-iknl 
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where [fc min ,fc max ] = 27r[/ min , / max ]/c is the wave num- 
ber, i.e., the frequency interval under consideration and 
£ is a geometrical length, and compared to the FT of the 
semiclassical trace formula, Eq. Q . From Eq. ^ follows 



— ikj n£ 



FT{p} , 



(10) 



and \p(£)\ — |FT{pfl uc (A:)}| is called the length spectrum. 
The experimental length spectrum \p(£)\ is computed by 
inserting the resonance frequencies fj and the widths 
r, obtained from the measured frequency spectrum into 



Eq. ([10]) . The FT of the contribution of a single PO to 
the trace formula [Eq. ^fj] is 



Pp,r(£) — o-p,r smc[nAk(r£ p — £)/2] 



(11) 



with sinc(x) = $m(x)/x and 

p ' r ~ 7r|det(M;-l)|i/2 K P e ' [iZ) 

where Afc = fc max - fc min and k = (fc max + fc min )/2. Fur- 
thermore, each PO is counted several times in Eq. ^ de- 
pending on its symmetry with respect to the two mirror 
symmetry axes of the stadium and with respect to time 
reversal. For each symmetry that is lacking for a certain 
PO, its contribution is counted twice [59Tf6~T] . which leads 
to a symmetry factor s p € {1, 2, 4, 8}. It should be noted 
that the FT of the DOS is different from that of the 
transmission amplitude Sba : while the former (i.e., the 
length spectrum) is connected to the POs of the billiard, 
the latter is proportional to the propagator from antenna 
a to antenna b and is related to the classical trajectories 
between the antennas 39. 621. 



IV. THE EXPERIMENTAL LENGTH SPECTRA 

Figures [4] and [5] show the length spectra for the stadia 
SI and S2, respectively. In each figure the top graph 
shows the experimental length spectrum. Both are quite 
similar on a qualitative level, and are therefore treated 
in parallel here. Both length spectra feature two major 
peaks, which are located at £ = 1.262 m and 1.285 m in 
the case of SI and and at £ = 1.083 m and 1.104 m in 
the case of S2. Most other peaks are not higher than the 
oscillating background, which is as large as \p(£)\ ~ 5 in 
both cases. The bottom graphs of Figs. [4] and [5j show 
the FT of the corresponding trace formula. The vertical 
lines in the bottom graphs indicate the lengths £ po and 
the modulus of the amplitudes, 



A,, 



(13) 



of the POs used in the calculation of the trace formulas. 
All POs with up to 12 reflections at the boundary were 
used. The POs are increasingly dense for lengths close to 
the circumference U, but their amplitudes drop fast for 
£po U because they become more and more unstable. 
Several families of these whispering gallery type POs with 
up to 50 reflections were considered additionally. A full, 
systematic analysis of these POs and their contributions 
[33] . however, was forgone since the contributions to the 
trace formula that are relevant in the present work arise 
from POs with lengths smaller than those of the whisper- 
ing gallery orbits. The overall shape of the experimental 
length spectra and the FTs of the trace formula agree 
well. The vertical lines in the bottom graphs indicate 
the lengths of the POs. They clearly show that the ma- 
jor peaks in the length spectra in fact result from the 
interfering contributions of several POs with similar am- 
plitude Ap^ r each. According to Eq. (11), the contribu- 
tion of a single PO has a width cx 1/Afc. With / m ; n = 5 
GHz and / max = 20 GHz here, 1/Afc = 3.2 mm, which 
is larger than the typical length difference of neighbor- 
ing POs in the considered length regime. A much larger 
frequency interval would be needed to resolve the POs, 
however, we are limited by $id i n the experiment. 
Though the shape of the experimental length spectra and 
the trace formula predictions agrees well for the major 
peaks, the peak amplitudes of the former are about four 
times smaller than those of the latter. For the most part, 
this can be attributed to the large number of missing 
resonances in the measured frequency spectra. A total 
number of 137 resonances in the range of 5-20 GHz was 
identified for stadium SI, and 138 resonances in the same 
frequency range for S2. The total number of resonances 
in a given frequency interval can be estimated from the 
Weyl formula [T7| 

N Wcyl (k) = ^-k 2 + P(n)^k + 0(l), (14) 

47T 47T 

where N^ ey i(k) is the number of resonances up to wave 
number k, A is the area, and U is the circumference of 
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FIG. 4. Length spectrum for the stadium SI. The top graph shows the experimental length spectrum, and the bottom graph 
the FT of the trace formula. Note the different scales of top and bottom graphs. The vertical line in the top graph indicates 
the path length h defined in Eq. (|3j) , the gray bar its error, and the arrow the circumference U of the stadium. The vertical 
lines in the bottom graph indicate the lengths and the amplitudes Ap lT of the POs [Eq. ( |13[ )] used for the calculation of the 
trace formula. Some POs are shown as insets. 
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FIG. 5. Length spectrum for the stadium S2. The top graph shows the experimental length spectrum, and the bottom graph 
the FT of the trace formula. Note the different scales of top and bottom graphs. The vertical line in the top graph indicates 
the path length li , the gray bar its error, and the arrow the circumference U of the stadium. The vertical lines in the bottom 
graph indicate the lengths and the amplitudes A p , r of the POs. Some POs are shown as insets. 



the resonator. The term f(n) accounts for the bound- 
ary conditions at the dielectric interface and is given by 
Eqs. (28) and (29) in Ref. [T7]- Its numerical value is 
f fa 1.02 for the Teflon stadia considered here. Accord- 
ing to Eq. ( |l4"| , 3507 resonances are expected to exist 
in the regime of 5-20 GHz for stadium SI and 2643 for 
S2, so only 3.9% and 5.2%, respectively, of all resonances 
were identified in the measured spectra. Therefore, in- 
vestigations of the statistical properties of the measured 
resonance frequencies and widths are practically impos- 



sible. The reason for the large number of unobservable 
resonances is their short lifetime due to the large radia- 
tion losses. A further possible reason for the discrepancy 
between the experimental and semiclassical peak ampli- 
tudes is that corrections to the Fresnel coefficients for 
reflections at curved interfaces are needed [63]. This is 
further investigated in Sec. El 

A closer inspection of Figs. [4] and [5] reveals that the POs 
contributing to the major peaks of the length spectra 
are all confined by total internal reflection (TIR), i.e., 
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the corresponding angles of incidence at the boundary 
are all larger than the critical angle for TIR, a cr ;t = 
arcsin (1/n) s» 45°. The trace formula also predicts peaks 
corresponding to POs not confined by TIR, e.g., close to 
1.20 m in Fig. [4] and close to 1.0 m in Fig.[5j but the 
experimental length spectra only feature an oscillating 
background and no real peaks in these length regimes. 
This is also the case for the BBOs: no significant peaks 
are seen at 0.6 m (primitive BBOs, not shown here) or 
1.2 m (first repetition). Summing up, the comparison 
between the experimental data and the trace formula 
for the two chaotic stadia yields similar results as in 
Ref. |21j . where regular dielectric microwave resonators 
were investigated: there is a good qualitative agreement 
between the experimental length spectra and the trace 
formula predictions, the peak amplitudes of the former 
are smaller than those of the latter due to a large number 
of missing resonances, and the experimental length spec- 
tra show no peaks corresponding to POs not confined by 
TIR. 

The vertical lines and the gray bars in the top graphs of 
Figs. [4] and [5] indicate the path length corresponding to 
the subset of equidistant resonances [Eq. ([3])] and its er- 
ror, respectively. For both stadia, I is close to the lengths 
of POs that are confined by TIR, which accounts for the 
relatively high quality factors of the equidistant resonan- 
ces. However, I cannot be identified with a single, specific 
PO, since there are several POs with similar amplitudes 
A Pj r within the error bars of I. Therefore, it cannot be 
clarified whether the set of equidistant resonances can be 
attributed to scarred states or not. In fact, in Ref. [ST] it 
was found that there is not necessarily a connection be- 
tween a family of equidistant resonances and a single PO, 
and in Refs. [TTMl3| . resonant states of dielectric stadia 
were investigated and scarred states localized on several 
POs each were found. So the sequences of equidistant 
resonances in the spectra of SI and S2 might also be lo- 
calized on a set of POs, but further investigations are 
necessary to decide this. 



V. COMPARISON WITH NUMERICAL 
CALCULATIONS 

Since only a small part of all resonances can be iden- 
tified from measured frequency spectra, we solved the 
Helmholtz equation [Eq. Q] for stadium SI numerically 
using a boundary element method (64] also used, e.g., in 
Refs. QH |5D]. In the range of 1.0-13.2 GHz a total of 
1648 eigenvalues kj was found. This is in good agree- 
ment with the number of states predicted by Weyl's law, 
Eq. ( 14 ) . Only 72 resonances were detected experimen- 



tally in the same frequency range, compared to the 137 
resonances found up to 20 GHz. A survey of the com- 
puted eigenvalues is shown in Fig. [6] where the imagi- 
nary parts of the computed kj are plotted with respect 
to the real parts as dots. The resonance frequencies and 
widths extracted from the measured frequency spectrum 
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FIG. 6. Real vs. imaginary parts of the eigenvalues for sta- 
dium SI in frequency units. The x marks are the experimen- 
tal data, and the dots the numerically calculated ones. 



are indicated as x marks. They cover only a narrow 
strip of the complex plane close to the real axis, i.e., only 
modes with relatively high quality factors were found 
experimentally. The calculated and measured eigenval- 
ues agree only roughly, which is attributed to the un- 
certainties in the determination of the resonance widths 
and the index of refraction. The length spectrum for 
the computed eigenvalues is shown in Fig. [7] as dashed 
line, that for the experimental length spectrum as solid 
line. They look different from those in Fig. [4] because 
a smaller frequency /wave number range Ak is consid- 
ered here, leading to broader peaks [cf. Eq. (fTT])]. The 
experimental length spectrum shows peaks with ampli- 
tudes about three times smaller than that of the numeri- 
cal one due to the small number of measured resonances. 
The dotted line is the FT of the trace formula for sta- 
dium SI. It shows good qualitative agreement with the 
numerical length spectrum, even though its peak ampli- 
tudes are somewhat larger. This deviation cannot be 
explained by missing resonances, since according to the 
Weyl formula only a very small number of resonances, if 
any, is missing in the set of numerically calculated reso- 
nances. A possible explanation is that the trace formula 
is not very accurate for POs with an incidence angle close 
to the critical one and that, thus, higher order correc- 
tions are needed for such POs (J7]. It is known that the 
Fresnel reflection coefficients at curved interfaces must 
be modified depending on the ratio of the radius of cur- 
vature R and the wavelength [B3], and it was suggested 
in Ref. [5U] to replace the ordinary Fresnel reflection co- 
efficients [Eq. Q] in the trace formula by ones with a 
curvature correction, 
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FIG. 7. Length spectrum for the stadium SI evaluated in the frequency range / = 1.0-13.2 GHz. Shown are the length 
spectrum of the numerically calculated eigenvalues (dashed line), the experimental one (solid line), the FT of the trace formula 
with the ordinary Fresnel reflection c oeffi cients [Eq. (JsJ)] inserted (dotted line), and that obtained with the Fresnel reflection 
coefficients for curved interfaces, Eq. (151 (dash-dotted line). The arrow indicates the circumference U. 
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FIG. 8. Modulus of the Fresnel reflection coefficient with re- 
spect to the angle of incidence 8 for n = 1.425. Shown is 
the modulus of the reflection coefficient incl uding a curva- 
ture correction, jr curv |, according to Eq. (151 for kR = 3.1 



corresponding to / = 1 GHz with R = 0.15 m (solid line), 
for kR = 6.3 (2 GHz, dashed line), for kR = 15.7 (5 GHz, 
dotted line), and for kR = 40.9 (13 GHz, dash-dotted line). 
The gray solid line shows the modulus of the ordinary Fresnel 
coefficient, jrj, defined in Eq. Q. 



tain the dash-dotted curve in Fig. [7j As expected, its 
peak amplitudes are smaller than those obtained for the 
trace formula without including curvature corrections. In 
the case of the peak near 1.20 m, the peak amplitude 
matches that of the numerical length spectrum, and in 
the case of the peak near 1.27 m, it is even a bit smaller. 
Furthermore, near 1.20 m the peak of the trace formula 
with curvature correction is shifted to the right with re- 
spect to that of the ordinary trace formula. This can be 
explained by the frequency dependence of the phase of 
the curvature corrected Fresnel coefficients [23]. In sum- 
mary, the introduction of curvature corrections leads to 
a better agreement between the calculated length spec- 
trum and the trace formula predictions, however, devia- 
tions remain. Since the major peaks in the length spec- 
trum result from the contributions of several POs each 
(cf. Fig. [4|) , a detailed analysis of the influence of the 
curvature correction on the contribution of a single PO 
is not possible here. This is necessary also because cor- 
rections to the trace formula apply for POs that are not 
well isolated [301 US]. Therefore, further investigations 
of higher-order corrections to the semiclassical trace for- 
mula are needed. 



VI. CONCLUSIONS 



where to = nkR sin 0j, Hm (z) is a Hankel function of the 

first kind of order to, and H'^ (z) is its derivative with 
respect to the argument. The modulus of 7Y curv with 
respect to the angle of incidence for different values 
of kR is compared to that of the ordinary Fresnel coef- 
ficient in Fig. H] The deviat ions between the corrected 
and the ordinary Fresnel coefficients is considerable, es- 
pecially close to the critical angle and for low frequencies, 
that is, small kR. The ordinary Fresnel coefficients are 
recovered in the limit kR — > oo. The curvature corrected 
Fresnel coefficients were inserted into the trace formula, 
Eq. f6|, and its FT was calculated numerically to ob- 



The frequency spectra of two passive 2D dielectric mi- 
crowave resonators with chaotic classical dynamics were 
measured and analyzed. The corresponding length spec- 
tra were compared to the predictions of a semiclassical 
trace formula proposed in Ref. [IT]. Good qualitative 
agreement between the experimental length spectra and 
the predictions of the trace formula was found. How- 
ever, the peak amplitudes of the experimental length 
spectra were systematically smaller than those predicted 
by the trace formula, which is mostly attributed to the 
large number of missing resonances in the measured spec- 
tra. Furthermore, no contributions from POs that are 
not confined by TIR were found in the experimental 



length spectra even if predicted by the trace formula. 
Apparently, the experimentally observable, long-lived re- 
sonances correspond to the best-confined POs of the cor- 
responding billiard, even though individual resonances 
cannot be associated with specific trajectories like for 
systems with regular classical dynamics, e.g., the circle 
or square [21 . The same observations concerning the 
amplitudes of and the POs contributing to the exper- 
imental length spectra were made in experiments with 
regular dielectric resonators in Ref. |21j . so they seem to 
be valid for all types of dielectric resonators. Further- 
more, the measured frequency spectra contained subsets 
of equidistant resonances and we investigated whether 
these resonant modes are localized on a particular PO. 
There are indications that they are instead localized on 
several POs each, but no clear conclusion could be drawn 
from the data. Spectra containing one or more sequences 
of equidistant resonances were also observed for other di- 
electric resonators with low index of refraction (n s» 1.3- 
1.5) and various shapes [23 ESI 13 SSI El EE| • Thus, 
this seems to be a generic phenomenon for dielectric res- 
onators with low n, however, the origin and whether it 
is a common feature is not known. Finally, we computed 
the resonances of one of the investigated resonators nu- 
merically to obtain a complete spectrum of eigenvalues. 



The corresponding length spectrum has been compared 
to the corresponding experimental one and to the trace 
formula prediction. It was concluded that the trace for- 
mula overestimates the peak amplitudes of the length 
spectrum. We applied a curvature correction to the Fres- 
nel reflection coefficients entering the trace formula and 
found better agreement between the numerical length 
spectrum and the trace formula prediction. However, 
some deviations remained, and it would be worthwhile 
to systematically investigate further corrections to the 
trace formula of Ref. [T7] for a more complete under- 
standing of the length spectra of dielectric resonators. 
Corrections to the Fresnel reflection coefficients for finite 
wavelengths are of general interest for the understanding 
of the ray-wave correspondence in dielectric resonators, 
see for example Refs. [rJ31 IBTI IB"8"] . 
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